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Abstract. The charged fermion mass matrices are invariant under U(l) 3 symmetry linked to 
the fermion number transformation. Under the condition that the definition of this symmetry in 
arbitrary weak basis does not depend upon Higgs parameters such as ratio of vacuum expectation 
values, a class of two Higgs doublet models (2HDM) can be identified in which tree level flavor 
changing neutral currents normally present in 2HDM are absent. However unlike the type I or 
type II Higgs doublet models, the charged Higgs couplings in these models contain additional 
flavor dependent CP violating phases. These phases can account for the recent hints of the 
beyond standard model CP violation in the Bd and B s mixing. In particular, there is a range 
of parameters in which new phases do not contribute to the K meson CP violation but give 
identical new physics phases in the Bd and B„ meson mixing. 



1. Introduction 

Cabibbo Kobayashi Maskawa (CKM) mechanism of CP violation has been established as the 
dominant source of CP violation by observations at B-factories and Tevatron. These experiments 
have also shown some hints of physics beyond standard model [1—5]. One of them is CP violating 
observable Sj/^j^s OI — B~d mixing which is measured in time dependent CP asymmetry of 

— > J/ipK s by BABAR experiment at the Stanford linear accelerator center (SLAC) and by 
Belle experiment at KEK [6,7]. Value of Sj^ Ks obtained from fit using e^, Am g , \V c b\, a, 7 as 
inputs shows deviations from experimentally determined value at about 2a [8]. Another hint 
for new physics comes from determination of CP violating phase (f) s of B® — B® mixing. cp s 
is determined from analysis of time dependent angular distribution of decay products in flavor 
tagged decay B® — > J/ipcj> by CDF and DO experiments at Fermilab Tevatron collider through 
decay chain B® — > J/ijj(f), J/i/j —> <p — > K + K~ [9, 10]. SM predicted value for <p s is much 

small compared to the experimentally determined value. UTfit group has combined all the 
available experimental constraints on B s mixing and performed a model independent analysis of 
NP contribution to B° s - B° s mixing. From this analysis UTfit group has concluded that phase 
4> s of B® — Bg mixing determined from their analysis deviates from the SM prediction by about 
3cr [1, 11]. A similar analysis performed by CKMfitter group shows that the deviation from SM 
prediction is about 2.5 a [12]. 

Inconsistencies between SM prediction and the experimental determination in B® — B® and 
B® — B® mixing discussed here can be resolved if there is some new physics beyond SM which 
can give extra contribution to above mentioned CP violating observables. Here we explore the 
possibility of explaining these deviations in two Higgs doublet model (2HDM). The most general 
2HDM generates large flavor changing neutral currents (FCNC) for moderate Higgs mass. To 
avoid FCNC an additional discrete symmetry is imposed which results in type-I or type II 



2HDM. In these models the charged Higgs couplings do not provide any additional source of 
CP violation. Therefore they cannot explain possibly large CP violating phases in the neutral 
B mixings if confirmed in future. Hence we need to go beyond the type-I and type-II 2HDM 
to explain new CP violating phases. One way is to restrict the structure of FCNC couplings 
using flavor symmetries rather than eliminating them altogether [13-17] . This leads to models 
with suppressed FCNC and additional phases whose phenomenological consequences have been 
studied in [18-20]. A different aproach based on shared flavour symmetry is given in [21]. An 
intersting possibility is to have models without tree level FCNC but containing additional phases 
in the charged Higgs couplings. Here we wish to discuss such models motivated by the studies 
of flavor symmetries of mass matrix in 2HDM. 

Flavor symmetries are often applied to restrict the structure of Yukawa couplings all of which 
cannot be directly determined from experiments. The relation between the structures of mass 
matrices and symmetries was studied in [22-25]. It was shown by Lam that a symmetry can 
always be found under which an arbitrary neutrino mass matrix M v remains invariant [25]. 
Grimus, Lavoura and Ludl [22] showed that any Hermitian mass matrix MjAlj obtained form 
a fermion mass matrix Mf always possesses a symmetry Gf = U(l) x £7(1) x [7(1) and the 
corresponding G for the mass matrix of the Majorana neutrinos is Z 2 x Z 2 x Z 2 . Here we 
generalize this study to the non hermitian mass matrices in 2HDM. 

In the next section we will describe general 2HDM and the structure of FCNC and charged 
higgs interaction. Study of the symmetries of mass matrices in 2HDM and their consequences 
on FCNC and charged higgs interaction will also be given. In sction 3. phenomenology of 
this model in neutral mesons mixing will be given. In section 4. we will obtain constraints on 
NP paramateres in present case from the neumerical analysis of B® — B® mixing and charge 
asymmetry of like sign dimuon events. Last section will present a summary. 

2. Mass matrix symmetries and 2HDM 

Yukawa couplings for 2HDM are given as 

-£y = Q' L {T ld (j>i + T 2d (p 2 )d' R + Q' L {Ti u (j)i + T 2u (j) 2 )u R + H.c. , (1) 

where Ti q (i = 1,2; q = u,d) are matrices in the generation space. (f>i t2 denote Higgs doublets 
and 4>i = iT2<t>i- Q\l re f er to three generations of doublet quarks and primed fields in the above 
equation refer to various quark fields in the weak basis. The neutral component of a specific 
linear combination of the Higgs fields (f> = cos (34>i + sin f3e~ l6 '<p 2 is responsible for the mass 
generation 

M q = v(cos pV lq + sin PT 2q e i(> «) = V qL D q Vj R , (2) 

where (0°) = wcos/3 ; = vsin(3e l6 with v ~ 174 GeV and 6d = — 9 U = 9. The matrices 

V q L,R diagonalize M q . We define the matrices S q L,R as 

S qL ,R = V qL , R P q V^ R , (3) 

It can be seen that the mass matrices have invariance given by following equation. 

S\ L M q S qR = M q . (4) 

Where P q = diag.(e iai % e ia2 " , e ias "). Here S qLtR define two different [7(1) x 17(1) x 17(1) 
symmetries G u and Gd for up and down quarks. We put a mild requirement on possible 
S q L,R namely that the form of S q L,R be independent of the parameters tan (3 and 9 which 
are determined entirely in the Higgs sector. With this requirement and from eq.(2) of mass 
matrices, we get 

si L r iq s qR = T iq i = l,2. (5) 



This shows that individual Yukawa couplings should also respect the symmetry. Let us 
parametrize Ti q as 

r iq = v qL r iq v] R . (6) 

Eqs. (3,5) then imply 

Pq^iqPq = ^iq ■ (7) 

If G u ,Gd refer to the full U(l) x U(l) x U(l) symmetry with totally independent ai q then the 
only non-trivial solution of eq.(7) is a diagonal Ti q for every i and q. Yukawa couplings are then 
given as 

r tq = v qLli y qR , (8) 

where ji q are diagonal matrices with complex entries. More general forms for Ti q are allowed 
if one demands invariance with respect to subgroups of G u x Gd- All the flavor violations are 
induced by Higgs combination (f)p = — sin/?0i + cos j3 4>2 e~ %e . The couplings of the neutral 
component 4> F are given as 

-t% = q L F q q R cf> F +Jl.c. (9) 

with 

F i = V^ L (- sin f3T lq + cos (3r 2q e id «)V qR 

= (-sin/37i g + cos/3 7 2ge^) , (10) 

where we have used eq.(8) to obtain the second line. It is seen that the FCNC matrix F q become 
diagonal along with the mass matrices and the tree level FCNC are absent. But the phases of 
{Fq)a cannot be removed in the process of making the quark masses real and remain as physical 
parameters. The charged component H + of (j>p correspond to the physical charged Higgs field 
and its couplings in our case are given by 

H+ 

-£>H+ = — (uiLVijiFdfjjdjR - Ui R Vij(F*)iidjL) + H.c. (11) 

The above couplings are similar to the charged Higgs couplings in 2HDM of type-I and II. In 
those models, (F q )u are proportional to the corresponding quark masses rm q and are real. Here 
(F q )a are general complex numbers which can provide new phases in the Bd jS — B^ s mixing. 

An interesting class of 2HDM without the tree level FCNC have been obtained from general 
2HDM by assuming that two Yukawa couplings T\ q and T2 q are proportional to each other [26]. 
In the present case, the two Yukawa coupling matrices are not proportional to each other but 
the tree level FCNC are still absent. The phases of {F q )u in the charged Higgs couplings are 
dependent on the flavour index i unlike in models of [26] which are characterized by universal 
phases one for the up and the other for the down quarks. If the diagonal matrices ji q and 72^ in 
eq.(8) are proportional to each other then the present class of models reduce to the one in [26]. 

If S u l 7^ SdL then neither the Yukawa interactions given in eq.(l) nor the charged current 
weak interactions remain invariant under symmetries of the mass matrices. This means that 
radiative corrections will not preserve [27] the structure implied by eq.(8). However as in case of 
the 2HDM of type-I and type-II as well as the aligned models of [26], the full Lagrangian of the 
present model is formally invariant under the fermion number transformation qn^R — >• e ta ^qiL,R 
accompanied by the change in the CKM matrix elements Vij — > e taiu Vije~ ta i d . As a consequence 
of this all the radiative corrections in the model would display structure similar to the one ob- 
tained in the Minimal Flavor Violating [28] models. 



3. Neutral Mesons mixings 

In this section we discuss phenomenological application of the model in new contribution to 
neutral meson mixing induced by the charged Higgs couplings in eq.(ll). The (F^u and (F u )a 
entering the H + couplings are determined by the diagonal Yukawa couplings 7^ which also 
determine corresponding quark masses, see eq.(8). We make a simplifying assumption that the 
first two generation quark masses and the corresponding (F q )u are small compared to the third 
generation masses and (F q )^. Hence eq.(ll) reduces to 

-C H+ « ^ (^3(^)336* " t R V 3j (F:) 33 d jL ) + H.C. (12) 

The charged Higgs contribution to the K° — K° mixing arise only from the second term. The 
phase in (-FIJ33 can be absorbed in the definition of H + . As a result the above Lagrangian 
does not generate any new CP violating phases in the K meson mixing as long as (F q )jj are 
neglected for j = 1,2. However it can lead to non-trivial phases in the B q — B q mixing since the 
charged Higgs exchanges in this case involve both (Fa) 33 and (F u )33 and their phases can not be 
simultaneously removed. More interestingly, the interaction in eq.(12) distinguishes between the 
d and s quarks only through the CKM factor and not through additional phases. This results in 
strong correlations among the CP violation in B s and system. This can be seen as follows. 

B® — B q (q = d, s) mixing amplitude can be parameterized in the presence of new physics 
contribution as 

(B q \H SM + U NP \B q ) = (B q \H SM \B q ) (1 + Kq e< P ) = \ (B q \H SM \B q ) \p q e~ 2i ^ + ^ , (13) 

where f3 q represent the relevant phase in case of the SM and (ft q are the charged Higgs induced 
phases. The box diagrams involving WH and HH lead to new physics contribution involving 
the charged Higgs H. The interaction in eq.(12) lead to the following effective Hamiltonian at 
the weak scale: 

H NP = C x q Ll n L q LliX b L + C 2 q- R b L q- R b L , (14) 

where Ci 2 are the Wilson coefficients. Taking matrix element of the above equation and 
comparing with the SM result leads to 



p np Att 2 



G 2 F M^ VB S (x t ) 



d - 5/24C 2 



M Bq \ 2 B 2q 



m b + m q J Bi q 



(15) 



where Gf,Mw respectively denote the Fermi coupling constant and the W boson mass. 
So(xt) ~ 2.3 for m t ~ 161 GeV. tjb ~ 0.55 refers to the QCD correction to the Wilson operator 
in the SM, -Big, 25 are the bag factors which enter the operator matrix elements in eq.(14) and 
MB q ,rrib,m q respectively denote the masses for the B q mesons for q = d, s, b quark and the 
d, s quarks. The Wilson coefficients Ci j2 are independent of the flavour q = d, s of the light 
quark in B q . Mild dependence of K q on q arise from the operator matrix element multiplying 
C*2 in eq.(15). This leads to two predictions: To a good approximation, (i) k& pa k s and (ii) 
(pNP ~ ^np 'p^jg i m pii es f rom eq.(13) that 



AM d _ AM| M 
A M s ~ AMP 



(16) 



where AM q denote the values of the mass difference between heavy and light mass eigenstates 
in B q system in presence of new physics. Equality of and k s as well as (j)^ p and (j)^ p along 
with eq.(13) also imply 

4>d ~ 4> s (17) 



The detailed phenomenological consequences of this prediction are already discussed in [29] in 
a model independent manner. It appears to be in the right direction for explaining the CP 
violating anomalies. In case of Bj, the unitarity triangle angle j3 

sin 2/3 = 0.84 ±0.09 

as determined [8] using the information from V c b,€K and is found to be higher than the 
value 

sin 2/3 = 0.681 ±0.025 

obtained from the mixing induced asymmetry in B — > J/ipKs decay. Since the latter measures 
/3 + (pd, the above information can be reconciled [29] with a negative 4>d ~ —10°. <pd ~ <p s then 
implies a sizable asymmetry [29] = sin2(/3 s — cj) s ) ~ 0.4 in B s decay as indicated by the 
analysis of UTfit [11] or the CKMfitter [12] group. 

In a recent study global fit to various obseravbles was carried out in different scenarios [30]. 
One of the scenario has the relation Kd = k s and 4>^ p = 4>^ p . The resulting fit in this scenario 
is better than the SM fit. Results from this fit are a comman NP phase 2<f)d = 2<f) s = —14.4^^2 
and sin 2/3 = 0.83 ± 0.05 at 2a. These results agree with our results obtained here. 



4. Numerical analysis 

In this section we discuss our numerical analysis in which we obtain constraints on model 
parameters from study of B® — B® mixing. New physics contribution to B® — B® mixing 
is parametrized by UTfit group as [11] 

CB 3 e-^ = ( 1 + (18) 



AAf s SM ) 

Also from eq.(13), C Ba = |1 + K s e~ 2i ^ P ~^^\ and <j> Bs = -|Arg[l + K s e~ 2i ^ P -^]. UTfit 
collaboration has obtained the allowed range of the parameters Cb s and 4>b s as [11] 

C Ba = [0.68, 1.51] , 

cj) Bs = [-30.5,-9.9] U [-77.8, -58.2] (19) 

We calculate Cb s and 4>b s in the present model neglecting first and second generation masses 
and corresponding F q u. The Charged higgs mass is allowed to vary in the range 100 — 500 
GeV while /3 and couplings F q ^ are varied in the allowed range. Values of Cb s and 4>b s which 
can be obtained in this case is shown in Figure(l). 




Figure 1. Values of NP parameter Cb s and (pB a which can be obtained in limit when first two 
generation masses and F q u , where q = u,d and i = 1,2, can be neglected 



We have also checked that our model is compatible with the recent determination of like- 
sign dimuon charge asymmetry of semileptonic b-hadron decays by DO collaboration [31]. The 
model discussed here can give extra contributation required to explain the deviation of the 
experimental value of like-sign dimuon charge asymmetry of semileptonic b-hadron decays from 
SM prediction. 

5. Summary 

Using the flavor symmetries which does not depend on the higgs parameter tan (3 and 6, 
we obtained a class of 2HDM in which FCNCs can be eliminated without imposing discrete 
symmetries. Unlike type - 1 and type-II 2HDM, charged higgs interaction in this model contains 
new phases which are flavor dependant. In the limit of vanishing first and second generation 
masses and corresponding couplings F q u where q = u,d and i = 1,2, K° — K° mixing does 
not get any new CP phase while the B® — B® and B® — B® mixing gets new phases. New 
contribution to CP violation in this model in B^ — B^ and B® — B® mixing are correlated and 
it is possible to explain CP violating anomalies in B^ and B s systems. The model parameters 
can also generate like sign dimuon charge asymmetry required by the experimental data. 
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